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O ■ Abstract 

> : 

^ ' The multiple scattering of high-energy particles in a thick target is formulated in an im- 



pact parameter representation. A formalism similar but not identical to that of Moliere is 
obtained. We show that calculations of particle beam broadening due to multiple Coulomb 
scattering alone can be given in closed form. The focus of this study is on whether or not the 
broadening of the Coulomb angular distribution prevents the retrieval of nuclear-interaction 
information from measuring the angular distributions of charged particles scattered from 
a thick target. For this purpose, we study multiple scatterings with both the nuclear and 
Coulomb interactions included and we do not make a small-angle expansion. Conditions 
for retrieving nuclear information from high-energy protons propagating through a block of 
material are obtained. 
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1 Introduction and basic formalism 

Understanding and calculation of the broadening of a particle beam when it prop- 
agates through a block of materials are important not only to make multiple-scattering 
corrections for cross-section measurements in physics experiments but also to many appli- 
cations such as radiography by means of high-energy protons. Many different theories of 
the multiple scattering of electrons by thick targets have been formulated in the past sixty 
years The theory of Moliere PP[7] has received extensive attention because it gives 

the best agreement with data concerning the broadening of Coulomb angular distribution, 
arising from the multiple scattering of charged particles from the atoms in thick targets. An 
excellent and succinct derivation of the Moliere theory was given by Mott and Massey [8] . 

As one deals with a large (almost astronomical) number of scatterers in a thick tar- 
get, the numerical aspect of the calculation becomes extremely demanding. One can obtain 
good results provided that careful approximations are carried out. Over the years, successful 
parametrizations of the broadening of angular distributions due to Coulomb multiple scatter- 
ing have been established [9]. For hadronic projectiles, such as protons, nuclear interactions 
also contribute to multiple scattering. However, this latter aspect has not yet received suffi- 
cient attention in the literature. In this work, our focus is, therefore, on effects of multiple 
scattering on angular distributions in the Coulomb-nuclear interference region and in the 
region where the nuclear interaction dominates. We formulate the multiple-scattering prob- 
lem in such a way that not only in the case with Coulomb multiple scattering alone can 
analytical evaluations of the beam broadening become possible, but also the calculations are 
greatly facilitated when both nuclear and Coulomb interactions are taken into consideration. 

For a very high energy (e.g., > 20 GeV) proton scattering from a single nucleus, 
the Coulomb cross section decreases rapidly with the increase of scattering angles in such 
a way that the Coulomb cross section is negligible with respect to the nuclear cross section 
already at scattering angles as small as several milliradians. However, the Coulomb peak is 
rapidly broadened by proton-atomic nucleus multiple scatterings. Clearly, for the purpose 
of extracting the forward amplitude of the basic hadron-nucleus strong interaction, one 
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should use targets as thin as possible and then employ Moliere-type theory to correct for 
the Coulomb contribution from the measurement [10] . On the other hand, in many practical 
applications the thickness of the "target" is often fixed by specific needs, which is by no means 
thin. It becomes, therefore, interesting to know how much nuclear information can still be 
learned from hadrons scattered from a thick target. Certainly, the feasibility of learning 
nuclear information can bring added value to probing materials with hadronic beams. In 
other words, will the broadening of the Coulomb angular distribution render impossible 
any study of the nuclear cross sections? One naturally expects that the survival of nuclear 
information, if any, depends on the target thickness, i.e., on the number of nuclear scatterers 
that a proton encounters in a block of material. We use our formalism to examine this 
question. 

In this section, after deriving the basic multiple-scattering formalism, we discuss the 
important bearing of an ability to retrieve nuclear information on many applications. In 
Section II, we show the broadening of angular distributions by Coulomb multiple scattering 
in closed form. The broadening of angular distributions by combined Coulomb and nuclear 
interactions are studied by means of semi-analytical models in regions of small momentum 
transfers (Section III) as well as large momentum transfers (Section IV). We find that it is 
possible to retrieve nuclear information from protons scattering from a thick target. Conclu- 
sions and suggestions are presented in Section V. 

It is well known that high-energy elastic scattering is basically forward peaked, 
which allows to a very good approximation to neglect the longitudinal momentum trans- 
fer. Glauber [11] has shown that, in an impact-parameter representation, every function 
0{q) of the tranverse momentum transfer q, whether O is an amplitude or a cross-section, 
can be parametrized in terms of a profile function 0(6) defined in the impact plane, 



where azimuthal symmetry is assumed, Jo is the Bessel function of the first kind, and q 
and b are the moduli of the transverse momentum transfer q and impact parameter b, re- 



oo 



0{q) = (27r)-2 / db exp{iq-b)(l){b) = {2n)-^ / bdb Jo{qb) (P{b), 
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spectively. Nonessential complications {e.g. spins, etc) are here understood. Conversely, the 
profile function results from the inverse Fourier transform and is given by 

oo 

(p{b) = J dqexp{-iq-b)0{q) = '^tt J q dq Jo{qb) O (q) . (2) 



Without loss of generality, we consider a target which consists of one kind of nuclei. 
For a thin target of thickness t and atomic density p the probability that a beam particle 
undergoes a scattering is = t p 0"*°*, where 0"*°* is the total cross section |12j. (The subscript 
1 denotes the single scattering.) The transmission probability is, therefore, given by p*''""* = 
1 — From the definition of the differential cross section (Ji{q), one obtains the sum rule 
^tot _ J dq 0--^(^qj_ The scattering probability density pi{q) is related to the differential cross 
section by |T2],[T3] pi{q) = t p ai{q). It is the probability that a particle experiences scattering 
in the direction q. The sum rule of leads to the sum rule = j dq pi{q). Most often, 
only the modulus q counts, o"i(g) = (Ti(g). Hence, a\"' = njd{q^)ai{q). 

In thick targets the beam can bounce forward from many nuclei and/or electronic 
clouds and/or different atoms. These multiple scatterings are incoherent because the scatter- 
ers are separated far apart with respect to the ranges of the screened Coulomb and nuclear 
interactions so that the scattering waves are already in the asymptotic region before the 
next collision occurs. Furthermore, the target is not crystalline on a macroscopic scale, and 
thus the distance between scatterers are largely random. One must, therefore, add proba- 
bilities (not amplitudes) coming from individual scatterings. If one splits the thick target 
with thickness T into a large number of thin targets each with thickness T/N, then 
the total multistep probability for a particle to be transmitted without any scattering is 

ptrans ^ ^ _ a{"^T p / N)^ . 

The differential probability density for just one scattering in this situation with many 
thin targets, each with thickness T/N., is 

PM = N aM {T/N) p (l - af'Tp/Nf'' . (3) 
This represents a combination of scattering from any one single layer and transmission 
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through the remaining {N — 1) layers without scattering. The factor N in front of the right 
hand side accounts for the N layers, obviously. 

For double scattering one must count pairs of layers and fold two single-scattering 
probability densities pi. Hence, 

N-2 



N 



(4) 



where 



<72{q) = j dq' ai{q') ai{q- q'). 



(5) 



For triple scattering, an identical argument induces the result 
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(6) 



where 



'^siq) = J dq' j dq" <Ji{q') Mq" ' q') Mq- q") 



(7) 



is a double convolution. Again, the first factor is the counting of all triplets of layers. 



It is useful at this stage to introduce the profiles 

oo 

0i(&) = j dqexp{-iq-b)ai{q) ^2tt j qdqJo{qb) ai{q), 



(8) 



with 01 (0) — crj°*, and, Vn, 

(f>n{b) ^ j dq exp(-ig • b) C7„(g) = 



(9) 



If we define $i as the profile of Pi, then we see that, 
Mb)^N(^')Mb) {l-ar'Tp/Nf-\ 



(10) 



and, more generally, the profile of P„ reads. 




(jf * T p 
N 



n N-n 



(11) 



The "total multistep profile" due to Pm = Y.n=i Pn can be written as 
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N 



N 



-Pn 



•i + -l^-'%<-^lf-(i-./iv) 
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(12) 



where we have introduced the parameter u = a\°^Tp. Like is dimensionless. However, 

it is not a probability for a thick target. Indeed, because the mean free path. A, of a beam 
particle is l/{pa\°^) , therefore, u = T/A represents the average number of collisions of a 
beam particle when it passes through a target of thickness T. Hence u can be a very large 
number. 



In the limit, N —>■ oo, we obtain. 



$m(&) = exp 



^tot 



— exp[—i') 



(13) 



Here, we emphasize the nonlinear action of u on $jv/(^)- The second term in Eg. (1131) comes 
from the limit of Pq, it also shows that u has the meaning of a beam decay rate in a target of 
thickness T. Note that this latter term is neglected in the Moliere theory [Ij as well as in its 
reformulation by Bethe [7] and by Mott and Massey [8]. Hence, we find a form similar but 
not identical to that of the Moliere theory for the final, multistep probability density Pm '■ 

Mb) 



Pm{q) = (27r) exp(— z/) dbexp{iq-b) < exp 



^tot 



(14) 



Again, we note that PmIq) depends nonlinearly on u. It is appropriate to underline the 
importance of the above-mentioned extra term in the present formulation. If T is small 
rather than large, T = t, then u is also small. The presence of this extra term reduces Pm 
to pi = tpai, as should be for thin targets. On the other hand, 0i(fe) — >■ as 6 — > oo. 



Consequently, the first term in the curly bracket has a limit, exp[ ] — 1, that exactly 
compensates the second term in the curly bracket, ensuring the convergence of the integration 
for Puiq)- 

It is interesting to note from Eq. (IT^ that 
Pif = J dqPM{q) = 1 - exp(-i.) = 1 - PMt , (15) 

where Pm* is the total transmission probability. This last equation is a sum rule for the 
multiple-scattering probability. Multiple-scattering differential cross sections can be related 
to probability densities by the general relation 

PM{q) = TpaM{q) , (16) 



in the same way cri(g) is to Pi{q)- In summary, three steps thus occur in this formalism: (i) 
Fourier transform thin target data o"i into their profile 0i, see Eq. ([H]); (ii) Find u = T pa\°^ 
and exponentiate Tp0i, see Eq. (1131) : (in) Fourier transform back into a probability 
distribution Pm, see Eq. ( !T4l) . 

It is reasonable to parametrize the single-scattering distribution as (Xi = o-ic + cin with 
(Tic being the Coulomb cross sections and (Ti„ the sum of cross sections of nuclear scattering 
and nuclear- Coulomb interference. The separation of ai into aic and ain induces the same 
for the profile function: 0i = (pic + (pin- The relation 



exp(/?0ic) exp{P(pin) - 1 = [exp(/50ic) - 1] + exp(/30ic) [exp(/30i„) - 1] 



(17) 



with (3 = iy/a{°* then leads to a split of Pm as the sum of two probability densities. 



Mc 



exp(— z/j 



db exp(zg ■ b) < exp 



^tot 



(1^ 



and 



Mn 



exp(— z/j 
(27r)2 



dh exp{iq ■ b) exp 



>lc 



(b) 



a 



tot 



exp 



u ■ 



'In 



ib) 



a 



tot 



- 1 



(19) 
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This allows a perturbative consideration of nuclear effects at those angles where Coulomb 
scattering dominates. Note that the exponent in Eq. (ITSl) contains a denominator 0"^°* and 
not cr^°*. Hence Pmc is proportional to a pure Coulomb process with an effective value of u, 
namely z/g// = (cr*°Vo"i°*) ^■ 



Alternately, at angles where nuclear scattering might dominate, the roles of (pic and 
>!„ can be interchanged to generate similar formulae, namely. 



P' 



exp(— i/j 



(2vr) 



dh exp{iq ■ b) I exp 



u ■ 



^tot 



1 



(20) 



and 



P' 



exp(— z/) 
(27r)2 



dh exp{iq ■ b) exp 



^tot 



exp 



ib) 



^tot 



(21) 



One can parametrize the screened Coulomb interactions as CTic{q) = J2a Ca {,<f + '^a)""- The 
powers a, screening momenta Kq, and normalizations Ca are mainly functions of the charge 
Z of each individual nucleus. But (7i„ will depend on both Z and the mass number A. The 
global normalization of aM will also depend on the target thickness or the parameter v. 
Hence, the theory is essentially driven by three parameters of a thick target, namely, Z, A 
and V. Experimental measurements of ctm might, conversely, permit a determination of such 
three parameters when the nuclear nature of the target is unknown a priori, as is most often 
the case for radiographic studies where a recovery of A, besides Z, would be precious. Success 
will occur, however, only if multiple scattering does not spoil the information carried by A. 
This question is the main concern of the following sections. 



2 Broadening of angular distributions by multiple scatterings 

The mean-square width of the distribution ai represents a useful observable for 
the broadening of the cross-section distribution ai and can be defined by an integral 

oo 

= {<V J dqq'aiiq) = 2 7r«)-i J dq gVi(g), (22) 
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if it converges. The use of Eq. ([8]) and elementary properties of the Fourier transform allows 
us to write 



^^ = -&oU^^^j 0^(0)' ^''^ 

where the operator between brackets ( ) comes from a two-dimensional Laplacian in cylin- 
drical coordinate space. The same procedure gives the mean-square width of ctm(^ as 

where $m(0) = (1 — e~'^). Assume, for the sake of the argument, that cii(g) is a Gaussian, 
a,{q) = Cexp[-qy{2K% (25) 



where C is a suitable normalization and the parameter 1/k is the interaction range. For 
instance, if one discusses screened- Coulomb interactions, then 1/k is of the scale of an atomic 
radius. Then one obtains a\°^ = 2nC , = and 



= 2TrCK' exp(-6'K72). 



(26) 



From Eqs. (fT3|) and (!24|) one further obtains 



(27) 



and 

(,V = ^^^^. (28) 

1 — exp(— z/) 

Since v is large in general, the denominator is ~ 1. Consequently, the multiple scattering 
has broadened the mean-square width by a factor as might be expected from a Brownian 
motion in the transverse-momentum space. The multiplication of {q^)i by v also occurs if 
we start from a "polynomial Gaussian distribution" (g/fi;)^" exp[— g^/(2K^)]. This growth 
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rate is very general and can be viewed as one more version of the central limit theorem. As 
additional evidence, one finds that if is of the functional form cri(g) = C(g^ + /t^)^" with 
an exponent n > 2, then {q^)i = jiji — 2) and again {q^)M = ^ ~ ^^^) — 

In the following, we illustrate the broadening of the cross-section distribution in the 
case of a screened Coulomb scattering. We fit cti at small angles by a few terms of the sum 

--(^)=E7-^%w, (29) 

m>2 '^mJ 

where m can be half-integers as well as integers, and Cm and Km are fitting parameters. It 
follows that 

C 

< = ^ E _ l^l^im^l) ■ (30) 



For definiteness, we take two terms with m = 5/2 and 4, namely, 

/ X _ ^5/2 Cj 

^^^^^^ - + + (^2 + ^2)4 ■ (^^) 



Hence, 



tot _ 2^C'5/2 -kC^ , . 



Dividing both sides by o"*°*, we obtain 

+ q 6 tot = "5/2 + ^4 . (^ddj 



3^5/2 "^Ic* 3^4(71° 



Eq. (1331) shows that both 05/2 and 04 are dimensionless numbers between and 1. An 
advantage of using Eq. fl29|) is that its Fourier transform gives the profile function 0ic in 
terms of analytical functions which can be easily analyzed, i.e.. 



2vrC5/2(l + /t5/2&) exp(-K5/26) 'nCAb'^K^{K^h) 



^^^^^^ = + 24.3 ^ (34) 



1 n 



where is the modified Bessel function of the second kind. One verifies easily that 0ic(O) 
2vrC5/2/(3KL2) + -nCJi-iHl) = a\°^. The "Coulomb" multistep profile then reads 



Mc = exp < z/ 



a5/2(l + «:5/2&)e-'^^/^' + ^i^3(«4&) - 1 



— e 



(35) 



When 6 — > CX3, it is easy to verify that Eq. fl34l) induces exponential decreases with ranges 
K^y^g ^^^d kJ^. Numerical integrals with such integrands converge well. The final integral for 
the Coulomb cross section then reads 

1 



2nTp 



dbbJoiqb) $mc(&)- 



(36) 



Let K be an average between the two momenta K5/2 and K4, which are both atomic 
scales. It is now convenient to scale momenta and lengths as g = nQ and b = B/n. The 
dimensionless Q will be a few units or a few tens, if one wants to describe scattering angles 
moderately larger than the Coulomb peak. We also expect that the values of B contributing 
to the integral 



(^Mc{Q) 



exp(— z/j 

2lTK,'^Tp 



dBBJoiQB) 



X < exp 



1^05/2 



(37) 



should be mainly between and several units. In atomic units, all parameters k, 1^^/2/1^ and 
Hii/ Hi are of order 1. It remains to estimate the dimensionless magnitudes of 1^05/2 and v a^l%. 
From Eq. fl33l) . 05/2 and 04 are moderate fractions of 1. It is thus the large number v that 
drives the integrand. 

It is also convenient to write 



ai(g) = /c-Vi(Q) , 01 (6) = h.-'' 4>^{B) 



(38) 



with cri(Q) and (f>i{B) being dimensionless. Because we work with systems of atomic scale, 
we further set k to be 1, meaning that our primary scale is "atomic" . In this scale, all lengths 

1 1 



and momenta will, respectively, be given in units of atomic radius and its inverse. 



To show the shrinking of profiles by multiple scatterings we plot 0ic and $mc in Fig- 1 
as the "crosses" and solid curves when aic = {Q"^ + and, respectively, as the "circles" 

and dashed curves when uic = {Q^ + As one can see, the dashed and solid curves do 
decay faster than their respective single scattering partners. The Fourier transform of $mc 




Fig. 1. Crosses: <^ic(S) for au{Q) = + Solid curve: the corresponding multistep $mc(-B) 

ii u = A. Circles and dashed curve: (t)ic{B) and ^Mc{B) (with v = A) for ctic = (Q^ + 1)"^. All 
profiles normalized to 1 at B = 0. 

then leads to the expected broadening of gmc, as shown in Fig. 2. (For graphical convenience 
we used v — Am. Figs. 1 and 2, which is much smaller than physical v but is demonstrative 
enough.) 
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Fig. 2. Crosses: aidQ) = (Q^ + l)'^/^. Solid curve: the corresponding (Tmc(Q) for = 4. Circles: 
(TiciQ) = {Q^ + Dashed curve: the corresponding aMc{Q) for = 4. All cross sections 

normalized to 1 at Q = 0. 

3 Loss of nucleeir information at small momentum transfers 



In this section, we introduce a semi-realistic model for ai{Q) which contains "nu- 
clear" information. We will investigate (a) changes of normalizations and not just shrinking 
or dilation of shapes of (j{Q), and (b) how nuclear information may become lost. We will, 
therefore, illustrate the blurring of signal through the study of various relevant quantities, 
such as Pm, Pmc-i Pum P'mci Pmu- analytical models to ensure that the blur- 

ring comes from physics and is not a result of numerical imprecision. A good analytical 
model must satisfy the following constraints: (i) positivity of the sum of and (Ti„; (ii) 
big contrasts between maxima and minima; (iii) analyticity in both the momentum and the 
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impact parameter representations; and (iv) significant differences between tlie atomic and 
tfie nuclear scales for profiles. 



Let 

<Jic(Q) = (1 + Q')-' (1 + gVlOO)-\ (39) 

This is qualitatively realistic, because the factor, (1 + <5^)~^, represents a screened Coulomb 
scattering. The additional, artificial factor, (1 + Q^/lOO)"^, is here just for the convergence 
of ((5^)ic- Then we further use, with (7i„(0) = 0, the following semi-realistic (7i„: 

CTiniQ) = 11 Q765060976287632087746874800000000 
X exp(l/2 - gVSO) [24684880296681586800 ei (3600 - 169 + 
-175731507577476 (Q^ - 144)^ (491485925 - 20593402 + 215573 Q*) 
+40919125 e W (Q^ - 25)^ (3973881778272 - 103796385841 + 463199137 Q^) ] (40) 

The quality of the model with respect to the requirements (ii) and (iv) are evidenced by Fig. 
3. 

As it is allowed by the split of ui into a "Coulomb" part (die) and a "nuclear" part 
{(Jin), our ain can be positive or negative, as long as ai remains positive. Our (Ti„ was fine 
tuned to create four clear "nuclear" signals, namely two maxima of ai near Q — 8 and 
18, and, as signatures of interferences, two sharp minima at Q — 5 and 12. Furthermore, 
we adjusted its parameters so that the maxima do not exceed ~ 1% of the forward peak 
of (Jic- Note also that our model ain has only two maxima and, thus, carries no "nuclear" 
information for Q > 40. This is designed to track whether or not the maxima, if they survive 
the blurring of angular distribution by multiple scatterings, would migrate towards larger 
values of Q. The log^o Ci and log^^Q aic of our toy model are shown as functions a{Q) in Fig. 
3. It is trivial to deduce (7i„ visually. 

The corresponding profiles read, in closed forms, 

, / X f-Ko(B) Ko(lOB) BKi{B)\ , , 

0icS)=2OO7r ^^ + ^^ ^ + , (41 

^^"^ ' \ 9801 9801 198 / ' ^ ' 
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Fig. 3. logio '^1 (solid curve) and logio (^ic (dashed curve) as functions of Q. 

and, with U = -6^2/25 ^ ^ _gll9/50^ 

(f)in{B) = 11/2081951241204226807899993600 tt exp(l/2 - 25B'^/2) x 
152587890625 (37882968282999233748 + 24684880296681586800 

+18953703386795125 v) - 4638671875 (1727605574138750181696 + 
1184874254240716166400 u + 973536146177791625 v) + 2734375 x 

(1318896762458059001549772 + 935772074343960364976400 w + 
815387774668679551625 t;) - 625 B^ (987420216355162578690908328 + 

716741544586342715769420000 M + 649720700222675824315625?;) + 
2800 (12734093638656544401340038 + 9458232286390484639445000 m 
+8603026609006638128125?;) ~ 8 (53363090971589265717153336 + 

40407914801652923512260000 u + 34558380311245783728125 v)\ . (42) 

These profile functions are shown in Fig. 4. The width of 0i„ is significantly smaller 
than that of (pic, as one should expect when comparing a "nuclear" profile to an "atomic" 
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one. A geometrical ratio of widths might be ~ 10~^ or even ~ 10~^, but the model ratio 
we choose, between ~ 1/5 and ~ 1/10, is sufficient for a pedagogical study and much more 
convenient numerically. 




Fig. 4. logxo '/'I (solid curve), log^o 0in (dotted curve), and log^Q epic (dashed curve) as functions of 
B. 

This choice of "data" gives, after a numerical implementation of Eq. ([Hj), the total 
multistep probability distributions shown in Fig. 5. The solid curve is the same as that 
in Fig. 3, namely, log^g ^i- The dashed, linked-crosses, and linked- circles curves represent 
logio Pm z/ = 4, 9 and 16, respectively. The result is striking, on two counts: (i) the 
forward peak is more and more damped, the distributions extending more and more towards 
larger momenta, and (ii) the "nuclear information", whether minima or maxima, becomes 
rapidly blurred beyond recognition. Furthermore, the broadening of distributions does not 
seem to push much residual information towards larger momenta. The broadening process 
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Fig. 5. The dependences on Q of logio '^i of single scattering (solid curve) and logio Pm of mul- 
tiple scattering when u = 4 (dashed curve), 9 (linked-crosses curve), and 16 (linked-circles curve), 
respectively. 

is also confirmed by the behavior of the component Pmc of Pm, shown in Fig. 6. In our 
model (T*"* ~ 5.7 and (T*°* ~ 3.0. We chose a large nuclear contribution, (tJ°* ~ 2.7, in order 
to emphasize nuclear effects. However, at Q < 40, even this exaggerated nuclear information 
did not survive multiple scatterings. 

In Fig. 7, we show the various probability distributions PuniQ)- We note again that 
multiple scatterings wash away nuclear information. A similar feature is also seen in the -P/vfn 
given in Fig. 8. Besides the damping and information loss which are evident from Figs. 7 and 
8, we may stress a feature of Fig. 8, namely the transformation of "negative cross sections" 
into positive ones after multiple scattering. In order to create interferences, it was necessary, 
at the stage of making a model for cri„, to create negative values interfering with Uic. As has 




Fig. 6. Solid curve: log^g o"ic(Q)- Dashed, linked-crosses and linked-circles curves: logio Pmc{Q) 
for v = 4:, 9 and 16, respectively. 

already been pointed out, this is allowed as long as ai remains positive; there is a degree of 
freedom in modeling o"i„. The solid curve in Fig. 8 shows log^^Q |cri„|. One sees four arches, the 
first and the tiny third ones meaning negative values. Such "negative" arches disappear in 
the dashed curves representing Pm^- This disappearance justifies the use of models where ain 
can be not everywhere positive as long as ain + <Jic is everywhere positive, as was discussed 
after Eq. 

An advantage of our use of special analytical forms for the cross sections, Eqs. (l39ll4Up . 
is that such forms induce analytical profiles, Eqs. fl41|42l) . which in turn allow analytical forms 
for the multistep profiles, Eqs. ([13]) and (fT8H2T|). Values of (Q^) 

can then be easily obtained 

from the use of Eq. ([23]). The rates of broadening as functions of u can also be readily 
calculated. Fig. 9 shows how, at values of u smaller by several orders of magnitudes than 
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Fig. 7. Dashed, linked-crosses and linked-circles curves: logio PMn{Q) for = 4, 9 and 16, respec- 
tively. 



those estimated from geometric cross sections, the square-widths (Q^) of PmtijPmc ah^eady 
increase Hnearly with u. We have also noted a similar behavior of the widths of Pl[^ and 
P' 

In summary, the signature of nuclear information (diffractive oscillations in the dif- 
ferential cross section) in the region of small momentum transfers is washed away by the 
broadening of the angular distribution. This happens even with our model that has exagger- 
ated nuclear cross sections. In the next section, we examine if there exist momentum-transfer 
regions where the multiple scattering of the proton does not completely blurr nuclear signals. 
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Fig. 8. Solid curve: log^o l^ijiCQ)! for a nuclear signal. Dashed, linked-crosses and linked-circles 
curves: probability distributions log^o PiiniQ) when = 4, 9 and 16, respectively. 

4 Conditions for observing nuclear information 

Let r = t>e the ratio of the half-width of the Coulomb peak to the mo- 

mentum transfer at which the first minimum due to nuclear diffraction is observed for a 
thin target. At high energies, both g^^^ ^-^d q^.^^-^ occur at very small angles. Consequently, 
^ ~ ^1/2 /^mmi with the 6''s being the respective scattering angles corresponding to q^^^ and 
g^.^j. From Eq. (!28|) . it is reasonable to expect a rule, (g^)Af — {q^)i — ^^f/a (~ '^'^^^Lni)' 
hence that there is a critical value Vcrit — f'^^, above which nuclear signals will be obliterated 
by the broadening of the Coulomb peak. In other words, nuclear signals can only be observed 

at g ^ qminl 

for U > Pcrit- 

In Fig. 10, we show the elastic scattering differential cross sections of protons scattered 
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Fig. 9. The dependences on log^o ^ of log^g (Q^) foi^ Pm{Q) (solid curve), Pmc{Q) (dashed curve) 
and PMniQ) (dotted curve). Note that the slopes ~ 1 when >> 10. 

from a thin ^^'^Pb target at ~20 GeV, which we have calculated by using the method of 
optical model of Ref . [H] with a screened Coulomb interaction. The calculated cross sections 
exhibit the main characteristics of high-energy proton-nucleus scattering, namely, a narrow 
forward Coulomb peak and the diffractive oscillations at larger angles. Here, the diffractive 
pattern constitutes the nuclear signal. One notes that the first diffractive minimum lies at 
about ~ 6 milliradians. We have noted from our calculation that 6*1/2 for the Coulomb peak 
is of order ~ 0.003 milliradians. It is therefore reasonable to assume that, for high-energy 
proton scattering from nuclei, r is of order ~ 10^"^ or less in general. At most one might 
consider r of order ~ 10~^. Accordingly, although the geometric size of a nucleus is typically 
~ 10^ — 10^ smaller than that of its atom, the range of "nuclear information profiles" at 
high scattering energies may be taken ~ 10^ to ~ 10^ smaller than the range of the screened 
atomic profile, and possibly much smaller. While the model used in the previous section 



where ~ 0.5 < r <~ 0.1 is pedagogically justified, an analysis with a smaller r is in order. 




Fig. 10. Differential cross sections of p— ^'^^Pb elastic scattering at 20 GeV. 

We first introduce a model in which r = 10~^; the profile function (f)i{B) is the sum 
of a "Coulomb" term, 

MB)^BK,(B) (43) 
and a "nuclear" term. 



4 X 10-6 
1 + exp[800(S- 1/100)] ■ 



= . , „„„ronn.o W.nnM " (44) 



The profile gives a bare Coulomb cross section of the form aic{Q) oc (1 + Q^)"^ and the 
Woods-Saxon profile 0i„ makes, in practice, a window with range r = 1/100 indeed. The 
coefficient 800 in its exponent creates a "smoothed" Heaviside function. Both profiles are 
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normalized so that 0"*°* = 1, and crJ°*/o"*°* = 4 x 10~^ except for a negligible factor 1 + e~^. 
This cross-section ratio is quite compatible with the suggested by Fig. 10. Hence, the set of 
parameters given in Eqs. and fl44p is more realistic than that used in the previous section. 
The result for various angular cross sections ctmIQ; u), compared with the single scattering 
cri(Q), is shown in Fig. 11. An inspection of the figure shows that the Coulomb peak damps 
and spreads and the nuclear signal fades when u increases. The solid curve, representing 
logiQ cTiiQ), and the dashed curve, representing log^^g (Tm{Q) for u = 2000, exhibit somewhat 
similar oscillations. The dotted curve, corresponding to log^g '^m for u = 10^, hardly oscillates 
any more, i.e., nuclear signals are completely washed out. This confirms a loss of nuclear 
signal at low and moderate momentum transfers when u approaches Ucru ~ f' '^ = 10*^. 
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Fig. 11. The dependences on Q of log^Q cri (solid curve) and of log^Q cja/ for z/ = 2000 (dashed 
curve) and u = 10000 (dotted curve). All with r = 10~^. 



The Ucrit ~ r ^ rule is also seen in the previous section, where the use of ~ 0.5 < 
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r <~ 0.1 induces the loss of nuclear signal as early as v >~ 10. To further verify this rule, 
we use the same but use instead r — 0.5 x 10~^ in 0i„, namely, 

4 X 10"'^ 

'^'"^^^ " l + exp[1600(i?- 1/200)] ' ^^^^ 

The results are shown in Fig. 12. As we can see, the observation of nuclear signals is much 
improved; in agreement with the rule that r = 1/200 elevates Vcrit to a higher value, ~ 4 x 10^. 




Fig. 12. The dependences on Q of logio (solid curve) and of F = logio '^M for v = 2000 (dashed 
curve) and v = 10000 (dotted curve). All with r = 0.5 x 10~^. 

As a last test of the Ucru ~ rule, we keep 0ic(-B) = BKi{B) and let 

2x10"^ 

'^'"^^ ~ l + exp[8000 (5- 1/1000)] ' ^ ' 
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which is an obvious r = 1/1000 case. The results are shown in Fig. 13, where it is clear that, 
as expected, Ucrit occurs between 10~^ and 10~^. 
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Fig. 13. The dependences on Q of log^o cri (solid curve) and of log^g ctm for v = 10^ (dashed curve) 
and u = 10^ (dotted curve). All with r = 10"^ . 

Because u — Tpal°* ~ Tpalf, the existence of a Ucru induces a critical target thick- 
ness Tcrit such that the retrieval of nuclear signal is possible for target thickness T sufficiently 
less than Tcru] namely. 



T <T, 



^crit 



crit — f„f • 



(47) 



For p-208pb elastic scattering at 20 GeV, al°^ ~ 6.7 x 10^ mb = 6.7 x 10"^^ cm^. The 
density d and the atomic mass number A of lead are 11.3 g/cm^ and 208, respectively. 
Hence, p = {d/A)HAvog = 3.27 x 10^^ cm~^, with the Avogadro number Havoq = 6.02 x 10^^ 
[1/mole]. Because from Fig. 10 it is likely that r ~ 10^^ and our analysis indicates that 
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Table 1 

Target thickness T corresponding to u = 10^ for 20-GeV protons. 



Atom 


Q [g/cm J 


L 


A 


A 


p [iU cm J 


T [cm] 


Pb 


11.3 


82 


1.000 


208 


3.27 


4.6 


Cu 


8.9 


29 


0.250 


64 


8.39 


7 


Al 


2.7 


13 


0.086 


27 


6.02 


29 


Mg 


1.74 


12 


0.077 


24 


4.37 


44 


Be 


1.85 


4 


0.018 


9 


12.38 


67 



Vctu ~ 10^, then Eq. (H7|) gives T^^t — 46 cm. Hence, a nuclear signal can be retrieved at 
V < 10^, which corresponds to T < 4.6 = 0.1 Tcru [cm]. In so far as Q = 1 corresponds to 
6 ~ 0.003 milliradians, the survivor oscillation seen for u = 10^ in Fig. 13 between Q ~ 6000 
and Q ~ 12000, compatible with the expected period ~ 27r/r, would demand experimental 
measurements at angles of order a few dozens of milliradians at most. 

The proton-nucleus Coulomb cross section a\°J is oc Z^Rl ~ Z^/^ where Z is the 
target charge and Re is the root-mean-square radius of electric charge distribution in an atom 
with Re ~ ao/Z^^^ and being the first Bohr radius [T5]. Hence, one can estimate cr*°*p^ 
for proton scattering from a given nucleus A at 20 GeV by using ~ A o"j;°*p_pf, with 

the scaling factor A = (82/Z)~^/^. Hence, for a same v one has T{pA) = T{pPb)ppb/{pA^)- 
Results for a sample of atomic nuclei at z/ = 10^ are given in Table I. Of course, the price 
one pays in studying the nuclear signals that survive the multiple-scattering broadening 
is that one has to measure the angular distribution with good energy resolution at large 
proton scattering angles. In the case of 20 GeV incoming protons, the angles are about tens 
of milliradians, where the magnitudes of the cross sections are quite small. However, such 
measurements should be feasible with the currently available technology. 
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5 Conclusion 

The main matliematical and pliysical statement proposed by tliis work about multiple 
scattering consists in folding probabilities rather than scattering amplitudes. This is justified 
by the incoherence which is expected between the different scatterers of a thick, non crys- 
talline target. Simultaneously, an eikonal approximation, justified by the very high energy 
of the beam, allows a familiar impact parameter representation with profiles. Furthermore, 
small-momentum expansions [Tj are not employed in the formulation. As a consequence of 
such initial statements, a Poisson process is found, leading to an elementary formalism of 
convolutions and exponentiations in a context of Fourier-Bessel transforms. 

This Poisson process is nothing but a random walk in transfer momentum space. 
The central limit theorem is at work and the details of nuclear oscillations and interferences 
between Coulomb and nuclear scattering are blurred very fast as soon as the parameter 
V = Tpcr\°^, a measure of the number of collisions, exceeds a critical value of order r~^. Here 
r is the ratio of the range of the nuclear profile to that of the atomic profile. 

Below this critical value of z/, and at moderate and large momentum transfers (at the 
cost of very small elastic cross sections in the latter case) our conclusion is that some nuclear 
information remains observable. Such information is contained in oscillations of the multistep 
angular cross section crAiiQ) with periods ~ 27r/r, oscillations that are similar to those of 
the Bessel function, Ji{rQ), which typically represents pure nuclear diffractive scattering. 
Our model analysis shows that the characteristic distances between successive cross-section 
maxima and minima in the angular distribution remain essentially unchanged while each of 
these oscillations dampens as u increases. 

From the point of view of retrieving nuclear signals from protons traveling through a 
thick target, for which no sufficient attention was given in the literature, our work is more 
of a general feasibility study rather than a specific numerical evaluation. We have made use 
of analytical and semi-analytical models to bring out the basic features of the underlying 
physics. We believe that the positive feasibility concluded from this study will sustain tests 
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in detailed numerical applications. 

One problem which has not been solved in the present work, however, is to find an 
estimate of the v dependence of such periods ~ 27r/r. Our numerical evidence suggests that 
the dependence is not very strong, despite all the causes for a broadening of the signal, but 
our models and calculations lack the precision needed to tabulate such periods into functions 
of V. This effort is under consideration for an extension of the present work. 

In summary, below t'crit, which is of order with r being the ratio of the range 
of nuclear profile to the range of atomic profile, the nuclear signals can be retrieved from 
proton scattering from a thick target of thickness T < Tcrit- We suggest a conservative upper 
bound, namely, T < 0.1 Tcrit for practical considerations. 
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